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SPECTRA OF INFINITE GRAPHS WITH TAILS 


L. GOLINSKII 


Abstract. We compute explicitly (modulo solutions of certain alge¬ 
braic equations) the spectra of infinite graphs obtained by attaching 
one or several infinite paths to some vertices of certain finite graphs. 
The main result concerns a canonical form of the adjacency matrix of 
such infinite graphs. A complete answer is given in the case when the 
nnmber of attached paths to each vertex is the same. 


Introduction and preliminaries 


0.1. Graph theory. We begin with rudiments of the graph theory. For the 
sake of simplicity we restrict ourselves with simple, connected, undirected, 
finite or infinite (countable) weighted graphs, although the main result holds 
for weighted multigraphs and graphs with loops as well. We will label the 
vertex set V(r) by positive integers N = {1, 2,...}, {v}v£V = oj < oo. 

The symbol i ^ j means that the vertices i and j are incident, i.e., {i,j} 
belongs to the edge set ^(r). A graph F is weighted if a positive number dij 
(weight) is assigned to each edge {i,j} G ^(b). In case dij = 1 for all i,j, 
the graph is unweighted. 

The degree (valency) of a vertex n G V(F) is a number 7 ( 1 ;) of edges 
emanating from v. A graph F is said to be locally finite, if 7 ( 1 ;) < 00 for all 
V G V(F), and uniformly locally finite, if supp 7 ( 1 ;) < 00 . 

The spectral graph theory deals with the study of spectra and spectral 
properties of certain matrices related to graphs (more precisely, operators 
generated by such matrices in the standard basis {efcjfcgN and acting in the 
corresponding Hilbert spaces C” or = t'^(N)). One of the most notable of 
them is the adjacency matrix A(F) 


( 0 . 1 ) 


A(F) = 


\0‘ij\\ij=l^ — 


dij, {i,j]e£{T)] 
0, otherwise. 


The corresponding adjacency operator will be denoted by the same symbol. 
It acts as 


(0.2) A{T)ek = '^ajkej, A; G N. 

j^k 


Clearly, A(F) is a symmetric, densely-defined linear operator, whose domain 
is the set of all finite linear combinations of the basis vectors. The operator 
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y4(r) is bounded and self-adjoint in as long as the graph F is uniformly 
locally finite. 

Whereas the spectral theory of finite graphs is very well established (see, 
e.g., [1, 5, 6, 7]), the corresponding theory for infinite graphs is in its infancy. 
We refer to [16, 17, 21] for the basics of this theory. In contrast to the 
general consideration in [17], our goal is to carry out a complete spectral 
analysis (canonical models for the adjacency operators and computation of 
the spectrum) for a class of infinite graphs which loosely speaking can be 
called “finite graphs with tails attached to them”. To make the notion 
precise, we define first an operation of coupling well known for finite graphs 
(see, e.g., [7, Theorem 2.12]). 

Definition 0.1. Let T^, /c = 1,2, be two weighted graphs with no common 
vertices, with the vertex sets and edge sets V(rfc) and T(rfc), respectively, 
and let Ufc G V(rfc). A weighted graph T = Ti -|- r2 will be called a coupling 
by means of the bridge {ui, U 2 } of weight d if 

(0.3) v(r) = v(ri)uv(r 2 ), £:(r) = £:(ri)u£:(r 2 )u{ui,u 2 }. 

So we join r2 to Ti by the new edge of weight d between V 2 and vi. 

If the graph Ti is finite, V^(ri) = {I,2,...,n}, and V^(r 2 ) = we 

can with no loss of generality put ui = n, U 2 = 1, so the adjacency matrix 
A(r) can be written as a block matrix 

0 ...‘ 

0 ... ' 

0 ... 

If r 2 = Poo({aj}), the one-sided weighted infinite path, we can view the 
coupling r = Ti -|- Poo({aj}) as a finite graph with the tail. This is exactly 
the class of graphs we will be dealing with in the paper. Each such graph 
has a finite number of essential ramification nodes (a vertex v is an essential 
ramification node if 'y{v) > 3, see [2]). 


(0.4) 


A(r) = 


AiVi 




h = 


0 0 

0 0 
d 0 



N + 2 

-• 


The spectral theory of infinite graphs with one or several rays attached 
to certain finite graphs was initiated in [12, 13, 14, 18] wherein a number 
of particular examples of unweighted (background) graphs is examined. We 
argue in the spirit of [3, 4, 20] and suggest a quite general canonical form 
for the adjacency matrix of such graphs and supplement to the list of the 
examples. As a matter of fact, the algorithm applies not only to adjacency 
matrices, but to both Laplacians on graphs of such type. We also find 
a canonical form of the adjacency operator for graphs obtained from an 



INFINITE GRAPHS 


3 


arbitrary finite graph by attaching an equal number of infinite rays to each 
of its vertices. 


0.2. Jacobi matrices. Under Jacobi matrices we mean here one-sided in¬ 
finite matrices of the form 


(0.5) 


J — J{{hj'\, {uj}) 


-bi 

ai 

ai 

b2 


02 


bj G M, aj > 0, 


which generate linear operators (called the Jacobi operators) on the Hilbert 
space £^(N). The matrix 


( 0 . 6 ) 


Jn := 


1 0 
0 1 
1 0 


0 

0 

1 


called a discrete Laplacian or a free Jacobi matrix, is of particular interest 
in the sequel. 

Given two Jacobi matrices k = 1,2, the matrix J 2 

is called a truncation of Ji (and Ji is an extension of J 2 ) if 




( 2 ) _ ( 1 ) 
0!j 


j gN, 


for some q gN. In other words, J 2 is obtained from Ji by deleting the first 
q rows and columns. If J 2 = Jo, Ji is said to be a Jacobi matrix of finite 
rank or an eventually free Jacobi matrix. 

For the class of Jacobi matrices of finite rank the complete spectral anal¬ 
ysis is available at the moment (see [8, 11]). A basic in perturbation theory 
object known as the perturbation determinant [9] plays a key role. Given 
bounded linear operators To and T on the Hilbert space such that T — To is 
a nuclear operator, the perturbation determinant is defined by 

(0.7) T(A;T,To) :=det(/+(T-To)T(A,To)), T(A,Tq) := (Tq - A)-^ 


is the resolvent of operator Tq, an analytic operator-function on the resolvent 
set p(To). 

The perturbation determinant is designed for the spectral analysis of the 
perturbed operator T, once the spectral analysis for Tq is available. In 
particular, the essential spectra of T and Tq agree, and the discrete spectrum 
of T is exactly the zero set of the analytic function L on p{Tq), at least if 
the latter is a domain, i.e., a connected, open set in the complex plane. 

In the simplest case rank(T — Tq) < 00 the perturbation determinant is 
the standard finite dimensional determinant. Indeed, now 
p p 

(T - To)h = Y,{h, Tfc) {T - To)R{X, To)h = ^(/r, R*{X, To)ipk) A, 

k=l k=l 

SO L can be computed by the formula 

(0.8) T(A;T,To) = det \\6ij + {R{X,To)fii,pj)\\l.^,. 
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Our particular concern is Tq = Jq, the free Jacobi matrix. Its resolvent 
matrix in the standard basis in is given by 

-| 1^_ji I 

(0.9) r(^z +-, Jo^ = \\rij{z)\\’^j^i, rij{z) = ^ ^ —’ zGlD), 

see, e.g., [10]. If T = J is a Jacobi matrix of finite rank p, we end up with 
the computation of the ordinary determinant (0.8) of order p. 

It is instructive for the further usage to compute two simplest perturbation 
determinants for rank(J — Jq) = 1 and 2. 


Example 0.2. Let 


J = J{{bj},{l}) : bj=0, j^q, 

so J — Jo = {•,eq) bqCq. Hence by (0.9) and (0.8) 

(0.10) L{z) := l(^z + = 1 + hqrqq{z) = l- bqZ ^ . 

Similarly, let 

J = J({0}, {ttj }) : Uj = 1, j q, 

so J - Jo = (•, eq) {aq - 1) Cq+i + {■, Cq+i) {ttq - 1) Cq. Hence by (0.9) and 
( 0 . 8 ) 


( 0 . 11 ) 


L{z) 


1 + (Oq 

(«<? - 

1 + ( 1 - 


- l)rq^q+l{z) 
1) rq+l,q+l{z) 


{aq - l)rqq{z) 

! + («<?- ^)rq+l,qiz) 


In the Jacobi matrices setting there is yet another way of computing 
perturbation determinants based on the so-called Jost solution and Jost 
function (see, e.g., [21, Section 3.7]). Consider the basic recurrence relation 
for the Jacobi matrix J 

(0.12) an-lVn-l + bnVn + Onl/n+l = (^ + “) 2/^’ Z GO, UG N, 
where we put ao = 1- Its solution = Un{z) is called the Jost solution if 
(0.13) lim z~^Un{z) = 1, z G B. 

n^oo 


The function u = uq m this case is called the Jost function. 

The Jost solution definitely exists for finite rank Jacobi matrices. Indeed, 
let 

Uq-l-l — aq-^-2 — • • • — 1, ^g+1 — bqJf-2 — • • • — 0. 

One can put Uk{z) = z^,k = q+1, q+‘2^, ■ ■ ■ and then determine Uq, Ug-i,..., uo 
consecutively from (0.12). So, 


(0.14) 


aqUq{z) = Z'', 

Oq^lOq Uq_i{z) = Uq Z^^^ - bqZ'^ J- z'^~^, Og := 1 - 


etc., and in general 
k 

Uq-k{z) = ^ 
j=-k 


k = 0,l,...,q, PqjGR, /3q-k = l. 
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In particular, for ^ = 1 

(0.15) ai u{z) = aiz^ — biz + 1, 

and for (7 = 2 

(0.16) 0102 ^( 2 ;) = a 2 z"^ — {h 2 + hia 2 ) z^ + {ai + a 2 + hib 2 ) z^ - {hi + h 2 ) z + l. 

The relation between the perturbation determinant and the dost function 
is given by 

CXD 

(0.17) u{z) = -Hz), 

i=i 

see, e.g., [ 10 ], and such recursive way of computing perturbation determi¬ 
nants is sometimes far easier than computing ordinary determinants ( 0 . 8 ), 
especially for large enough ranks of perturbation. 


Example 0.3. Let J = J{{bj},{aj'\) be a Jacobi matrix such that 
bj ~ 0 ) j 1 ) ~ 1 ) j Q- 

We have Uq^j{z) = , j = 1, 2 ..., 

(^q^qiz) — 2^^, CLqUq—\{^z) — Otq z'^ J- Z^ , 

aqUq-2 = Oiq ( 2 ;^+^ -k z"^) + 2!^“^, 
and, by induction, 

2k 2 

(0.18) aqUq-k{z) = aq , k = l,2,...,q-l. 

Next, for (7 = 1 we have exactly (0.15), so let q > 2. The recurrence 
relation ( 0 . 12 ) with n = 1 gives 

aqu{z) + biaqUl{z) + aqU 2 {z) = (^ + “) 0.qUl{z), 
and so we come to the following expression for the dost function 
(0.19) aqu{z){z^ — 1) = aq {z — bl) — blQqZ^ + Qqz'^ + biz — 1. 
Similarly, for the Jacobi matrix J = J{{bj}.,{aj'\) with 
bj = 0 , j ^ q', ~ 1 ) d 7 ^ 1 

one has 


. . / N , 2:^'?+^ J- aiz‘^'^ ^ — ai2;^ — 2; 2 

(0.20) aitt(2;) = -bq -^---h ai2;^ J- 1. 

2 ;^ — 1 


For the Jacobi matrix J = J({0},{aj}) with aj = 1, j / 1, q, the dost 
function is given by 


^ + ai2;^'^ — q;i2;^ — 2 :^ 2 

(0.21) aiaqu{z) = Oq -- - -h ai2: J-1. 


The spectral theorem for finite rank Jacobi matrices provides a complete 
description of the spectral measure of J [ 8 ]. 

Theorem DS. Let J = J({ 6 j}, {aj}) be a Jacobi matrix of finite rank 

a^+l — ®g+2 — • • • — 1) ^g+1 — ^q+2 — • • • — 0, 

and u = tto(J) be its dost function. Then 
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• n is a real polynomial of degree degw < 2q (the dost polynomial), 
deg u = 2(7 if and only if Og / 1 . 

• All roots of u in the unit disk B are real and simple, u(0) 7 ^ 0 . A 
number Xj is an eigenvalue of J if and only if 

( 0 . 22 ) Xj = Zj + Zj £ (—1,1), u(zj) = 0 . 

• The spectral measure cr( J) is of the form 

N 

(0.23) a(J, dx) = (JaciJ, dx) + crd(J, dx) = w(x) dx + aj5(Xj), 

i=i 


where 

2,\u(e«W ' ^ = u’(z,)u(l/z,) ■ 

Note that |u(e**)p = Q(x), x = 2cost, Q is a real polynomial of the same 
degree as the dost polynomial u. 


The dacobi matrices arise in the spectral graph theory thanks to the 
relation for the adjacency matrix A(Poo({aj})) of the weighted path 

(0.24) A(Poo({a,})):= J({0},{a,}). 

In case of the unweighted path we have 

(0.25) A(Poo) = Jo. 

Note that the adjacency matrix for the finite (unweighted) path P^ with m 
vertices is the finite dacobi matrix J({0}, {!}) of order m. The spectrum of 
this matrix is well known [5, p. 9] 

(0.26) o'(Pm) = |2cos I 

I m + 1J i=i 

It follows from (0.4) that for an arbitrary finite weighted graph G 
(0.27) A(G + Poo({a,})) = 

Sometimes two-sided dacobi matrices 


AiG) Id 


(0.28) J = J{{hj]j(z'z,{aj]j(z'E) = 


a -I ho ao 

ao hi 
ai 


ai 

62 0-2 


bj G M, ttj > 0, which generate linear operators on the two-sided arise 

in the canonical models for certain infinite graphs. We say that such matrix 
has a finite rank if there are integers N^, N~ < such that 

bn = 11, ni[N-,N+], an = l, ni[N-,N+-l]. 

In this case one has a pair of dost solutions {u^} of the basic recurrence 
relation (0.12) (with n G Z) so that 

(0.29) Un(z) = z^, n > u~(z) = z~'^, n < N~. 
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The Wronskian of two solutions {fniz)}nez and {gn{z)}nez of (0.12) is 
defined as 


(0.30) [f,g] := an{fniz)gn+i{z) - fn+i{z)gn{z)) 

(the right hand side does not actually depend on n). We denote by 

(0.31) w{z) := [u'^{z),u~{z)], z G D 

the Wronskian of two dost solutions of (0.12). 

It is well known that the spectrum of a two-sided Jacobi matrix J of the 
finite rank is 


a{J) = [-2, 2] U ad{J), ^J) = 

and for the eigenvalues Xj the following relation (see, e.g., [22, Theorem 10.4]) 

(0.32) 3hefiZ):Jh=(^C + ^'jh^w{C) = 0, c e 

holds. So the eigenvalues of J are exactly the Zhukovsky images of zeros of 
wronskian w (0.31) in the unit disk. 

Example 0.4. Let 

(0.33) J = J{{0}j^i,{aj}j^z), aj = 1, j ^ 0. 

We now have N~ = 0, N~^ = 1. A simple computation gives 


w(z) = - QqZ, 

aoz 


and so 





0-0 > 1 , 


and (Td{ J) is empty for 0 < uq < 1. 

Example 0.5. Let 

(0.34) J = J({0}j6z, {aj}jez), aj = 1, j / ±1. 

We now have N~ = — 1, N~^ = 2. A simple computation gives 

, , o / «! 0—1 1 \ 1 

w{z) = aia-iz + -1-^ H-, 

Va_i ai aia_i/ a^a-iz 




and so to determine the discrete spectrum ad(J) one has to analyze the 
roots of the quadratic equation 

(0.35) Q(y) = Ay^-By+1 = 0, A := (a^j^-l)(af-l), B = af+a^j^-1. 


Precisely, each root yo of this equation in (—1,1) generates two symmetric 
eigenvalues 


A± 


:= ± 
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0.3. Canonical form for adjacency matrices. It is well known that each 
bounded and self-adjoint linear operator on a separable Hilbert space is 
unitarily equivalent to an infinite orthogonal sum of Jacobi operators. We 
suggest here a “canonical” form for certain infinite adjacency matrices and 
the algorithm of their reducing to such form. 


Theorem 0.6. Let A he a bloek matrix in ^, 

I'O 0 0 

U U 


(0.36) 


A = 


T* T 


Id = 


000 ... 
d 0 0 ... 


whereA= is a real symmetric matrix of order N, J = J({/3j}, {aj}) 

is a Jacobi matrix. Then A can he reduced to the block diagonal form 


(0.37) 



J 


where A is a real symmetric matrix of order n — k, for some 1 < k < n, and 
the Jacobi matrix J is the extension of J. In other words, there is a unitary 
operator U in such that 

(0.38) U-^AU = A®J. 


Corollary 0.7. Given a finite weighted graph G, the adjacency operator of 
the coupling T = G-|-Poo({aj}) is unitarily equivalent to the orthogonal sum 

(0.39) [/■M(r)c/ = F(r)© j(r) 

of a finite dimensional operator F(T) and a Jacobi operator J{T) which is 
the extension of J{{0},{aj}). IfF is unweighted, J{T) is of finite rank. 

We call F(r) a finite-dimensional component of the coupling F and J(r) 
its Jacobi component. 

In particular, the spectrum of F is 

(0.40) fT(F) = [-2, 2] U fTrf(F), fTrf(F) = w < oo, 

is a union of eigenvalues of F(F) and J(F). By Theorem DS, the spectral 
analysis of the graphs with tails amounts thereby to finding their finite¬ 
dimensional and Jacobi components, computing the Jost polynomial for 
J(F) and solving two algebraic equations, the characteristic equation for 
F(F) and the Jost equation for J(F). 


Remark 0.8. Given a finite graph G, one can attach p > 1 copies of the 
infinite path Pqo to some vertex v G V(G). Although the graph F thus ob¬ 
tained is not exactly the coupling in the sense of Definition 0.1, its adjacency 
operator acts similarly to one for the coupling. Indeed, it is not hard to see 
that 


(0.41) 


A(F) 


r^(G) 


0(0 Jo). 

2 = 1 


Hence Theorem 0.6 applies, and the spectral analysis of such graph can be 
accomplished. 



INFINITE GRAPHS 


9 


The algebraic equations which encounter later on cannot in general be 
solved explicitly. We can only determine how many roots (if any) they have 
in (—1,1), by means of the following well-known result (see, e.g., [19, p. 41]). 

Theorem (Descarte’s rule). Let a{x) = aox'^ an be a real 

polynomial. Denote by iJ,{a) the number of its positive roots, and p(a) the 
number of the sign changes in the sequence {ao,... ,an} of its coefficients 
(the zero coefficients are not taken into account). Then p(a) — /u(a) is a 
nonnegative even number. 

Acknowledgement. I thank Professor Yu. Samoilenko for drawing my 
attention to the spectral graph theory and putting forward some particular 
problems related to infinite graphs with tails and their spectra. 

1. Proof of the main result 

Before starting the proof we note that the algorithm of constructing a 
new basis suggested below is even more important than the result itself. To 
determine the spectra of graphs in the next Section we will have to carry it 
out by hand. The result of Theorem 0.6 only guarantees that the procedure 
can be accomplished with enough paper and patience. 

The algorithm applies in a number of situations beyond graphs with one 
tail (for instance, for the chain of cycles, the ladder with missing rungs etc.). 
We will elaborate on this topic in our forthcoming papers. 

Proof of Theorem 0.6. 

We construct an orthogonal basis {efc}fc>i in so that the matrix of A in 
this basis has the canonical form (0.37). We take e^ = e^, k = n,n + 1,..., 
so 

Aen-\-l — den T fdien+l T OlCn-|_2, 

ACn+k — ^k—l^n+k—1 T Pk^^n+k T CkfcCn+fc+lj k — 2, 3, ... . 

We want to find an orthogonal matrix B = (more precisely, the 

matrix with orthogonal columns) such that 

n—1 

(1.2) ek — ^ ^ bikOi, k — 1,2,..., 71 1. 

i=l 

The following notation will be convenient throughout the proof 

n—1 

(1.3) Wjk ■='^bikaij, j = 1,2,... ,77, A: = 1,2,... ,77 - 1. 

i=l 

1. Put 

n—1 

&i,n—1 •— i — 1, 2, . . . , 71 1, ^n—l •— ^ ^ 

i=l 

and so 

(1.4) — ^n—1 “1“ ^nn^n 4“ ^^n+1* 
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Next, 


n—1 


n—1 n 


^^n—1 — ^ ^ ^ ^ ^ ^ ^ ^ '^j,n—l^j 


2=1 

n—1 


2=1 j=l 
n—1 


i=i 


i=i 


■ ^n- 


i=l 


Note that Yli ~ Xli ^in-i > 0, for otherwise the original matrix already 
has the canonical form. Define 

(1.5) yi :=—-, bj^n-2 ■= Wj^n-l - yibj,n-l- 


ipn-i ,2 

2^j=l 


Then 


n—1 


Asn-l = en -2 + yien -1 + ZiSn, Zi := ^ b^ 


n—1 ’ 


2=1 


and orthogonality of the columns 

n—1 

(1-6) yy bi^n-2h,n-l = 0 

2=1 

follows directly from the choice of yi (1.5). 

If en -2 = 0, that is, bi^n -2 = 0 for all i = 1,... ,n — 1, the algorithm 
terminates, since 

ACn—l — 2/lCn—1 T ZiCm 

the subspace := span{ej}j>n-i is invariant for A and so is its orthog¬ 
onal complement C^_i = © Cn-i, which is finite dimensional. Thus, we 

can extend the basis {ej}j>n-i in ^n-i in an arbitrary way to the basis in 
the whole The latter is exactly the basis we are looking for. 

If 6^-2 / 0 we proceed to the next step. 

2. We have 

n—1 n—1 n n 

^^n—2 — ^ ^ &2,n—2-^^2 — ^ ^ &2,n—2 ^ ^ ^ ^ '^j^n—2^j 

2=1 2=1 j = l j = l 

n—1 n—1 

— ^ ^ '^j^n—l^j “1“ '^n,n—2^n — ^ ^ '^j,n—l^j 

j=l j=l 


iWn,n-2 = 0 by (1.6)). Put 
(1.7) 


X]i = l '^j,n-2bj,n-2 


2/2 := 


y^n -1 ,2 
^ 1=1 ©,n -2 


, 2^2 := 


Wj^n-2bj^n-l 


y^n- 1,2 
^i=l ©>-l 

and define 

(1.8) bj^n-3 ■= yJj,n-2- y2bj,n-2- Z2bj^n-l, j = 1, 2, . . . , n - 1, 

which leads to 

(1-9) Aen -2 = en-3 + y 2 en -2 + 2 : 2 en-i- 
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The orthogonality relations 

n —1 n —1 

(1-10) ^ ^ 3^^,n—2 — &z,n—3&i,n—1 — 0 

i=\ i=l 

stem from the definition of 2 / 2;^2 and (1-6)• 

Again, if 6^-3 = 0, we are done. If not, we proceed to the next step. 
3. As above we have in view of (1.14) 


Put 
( 1 . 11 ) 
and define 


n —1 


Aejj_3 — Wj^ji— 3 ^j — ^ 

i=i i=i 


E n —1 7 7 

1=1 3Pj,n—3 /_77 = 1 ^j,n—3^j,n—2 


n —1 


2/3 := 


Y^n-l72 

^1=1 ^j,n-3 


, 2:3 := 


l2 

A^i=l ^j,n-2 


( 1 ' 12 ) jj _4 .— Wj^YL—2 y2bj^n—2 ^2bj^n—l^ j — 1 , 2 , . . . , 71 1 , 


SO 

(1-13) Aen_3 = en-4 + 2/3en-3 + Z3en-2- 

As above, the orthogonality relations 

n—1 n—1 

(1.14) ^ ^ ^2,?!—4^2,n—3 — ^ ^ &2,n—4&2,n—2 — 0 

2=1 2=1 

stem from the definition of 1 / 3 ,Z 3 and (1.14). As for the last one we see from 
(1.12) and (1.14) that 

n—1 n—1 n—1 

^ ^ ^ 2 , n — 4 ^ 2 , n —1 — ^ ^ ^ 2 , n — 3 ^ 2 , n —1 — ^ ^ 

2 = 1 2=1 2 = 1 

n—1 /n—1 \ n—1 

— ^ ^ ^j,n—3 I ^ ^ bi^fi—l^ij I — ^ ^ ^j,n—3^j,n—1 — 0 
j=l \i=l / j=l 

since wj^n-i = &i,n -2 + yibj,n-i (see (1.5)). 

k. As soon as the process has not terminated, the pairwise orthogonal 
columns {bi,n-i}^Zi i ■ ■ ■ {bi,n-k}'iZi have been constructed, 

n—1 

(1.15) ^ ^ bi^Yi—pbi^n—q — 0, P Qj P'iQ — 1, 2, . . . , fe, 

2=1 

and the equalities 

(1.16) 

& 2 ,n—1^2,n—m+l“l“-2^2n—l& 2 ,n—m+2 — ^ 2 ,n—m+1? — 2, . . . ,/c, Z\ — 0 

hold. Next, we have 

n n—1 

■^^n—k — ^ ^ j^n—k^j — ^ ^ '^j^n-k^j ; 

j=l j=l 



3(Xji I 6i,n—1 
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as Wn,n-k = 0 ill view of (1.15). Again put 


(1-17) vk ■■= 


Xyj=l ^j,n—kbj,n—k Xjj = 1 ^j,n—kbj^n—k+1 


y^n-1,2 
L^j=l j,n—k 


y^n-1,2 
iLjj=l ^j,n—k+l 


and define 


(1.18) bj n—k—l ■ '^j,n—k Uk^j^n—k ^k^j^n—k+lt 3 1 , 2, . . . , 71 1 , 

SO 

(1.19) Ae-fi—k — Gfi—k—l “1“ Vk^^n—k “1“ fc+l- 

It remains only to verify the orthogonality conditions with the last column 
{hi^n-k-i}■ The first two 

n—1 n—l 

^ ^ bi^n—k—lbi^n—k — ^ ^ ^i,n—fc—l&i,n—fc+l — 0 


follow directly from the choice of yk,Zk- For m = l,2,...,A: — 2 


n—l / n—l 


^ ^ bi^n—k—lbi,n—m — ^ ^ U)i,n—kbi,n—m — ^ ^ I ^ ^ bj ^-k^ji 1 

i=l i=l i=l \i=l / 


as claimed. 



If Cn—fc —1 — 0 then Acfi-k — Uk^^n-kA ^fcCn—fc+i) ^n—k •— span{ej}j>n,—fc is 
A-invariant, and so is its orthogonal complement which is finite dimensional. 
Otherwise the algorithm can be extended to yet another step. 

When the algorithm continues till ei, we come to the orthogonal set 
{®i}j=i ia span{ej}i<j<n-i, so necessarily eo = 0. In this case the finite¬ 
dimensional component is missing. The operator A in the normalized basis 
{efc}fc>i, Cfc = ||efc||“^ e-k is given by the Jacobi matrix which agrees with the 
original matrix J from some point on. The proof is complete. 


2. Spectra of certain trees with tails 

As we have already mentioned in Introduction the adjacency matrix A(r) 
for the coupling T = G J- Poo({aj}) is (0.36). The algorithm suggested 
in Theorem 0.6 provides a way to implement the spectral analysis of the 
adjacency operator A(r) which consists of two stages. First, one has to 
apply the above algorithm by hand to obtain the canonical form (0.37) and 
to determine the finite-dimensional and Jacobi components of F. Second, 
the solutions of two algebraic equations (the characteristic equation for T(F) 
and the Jost equation for J(F)) provide the spectrum of F. 

Example 2.1. “A weighted star”. 

We begin with the coupling F = Sn{w) + where Sn{w) is a simple 
weighted star graph of order n + 1, n > 2, with vertices 1 ,..., ri of degree 1 
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and the weight of the edge {k, n + 1) is Wk, 1 < k < n. The canonical basis 
{efcjfcgN looks as follows. We put 

Cj := Cj, j > n + 1 => ^(r) Gj = Gj-i + Gj+i, j > n + 2. 

Next, let w := (wi,W 2 , ■ ■ ■, Wn), ||in|| = + ... + w'^, and let 


^n. • — 


W\ 


E 

i=i 


WkGj 


Then 


A{T)en+i = \\w\\en + en+2, A{T)en = ||ic||en+i. 

So the Jacobi subspace and Jacobi component of T are 
(2.1) £j = span{e,}j>„, J(r) = J({0}, {||u;||, 1,1,...}). 

To find the finite-dimensional component we construct an orthonormal 
basis in C” by means of a unitary matrix ^ with the specified 

last column 

n 

( 2 - 2 ) fj-='^ikjGk{n), ^kn 


k=l 


Wk 

Ikll 


k = 1,... ,n, 


where {efc(n)}^^^ is the standard basis in C"". Put 

(2.3) Gj := {fj,0,0,..j = l,...,n. 

The orthogonality relations {Gk,Gn) =0, l<A:<n — 1, give 

n 

(2.4) yl(r) Cfc = ^ ikjikn ■ Gn+l = 0. 

i=i 

Hence the finite-dimensional component F{T) = 0,i_i on the subspace 
spanjcj ■ So the canonical form is 

(2.5) yl(r) ~ On-i 0j({O},{||rc||,1,1,...}). 

The Jost polynomial is now given by (0.15) 

||rt;|| u(z) = (1 — ||rc|p) 2 :^ J- 1. 

Clearly, u > 0 for ||rt;|| < 1, and it has zeros inside (—1,1) if and only if 
11 re 11 > \/2- In this case the discrete spectrum is 

1 


(2.6) i7rf(5n(in)+Poo) = {0 (”-i),±(v1HF^ + 


-1 


)}■ 


For the unweighted star Sn the discrete spectrum is 

1 


(2.7) 


crrf(5n-kPoo) = jo^"' ± - 1 -k ^ I, n > 3, 


(Td{S2+F^) = {0^^'>}. 


This case is studied in [14]. 

Note that Sn is a complete bipartite graph, Sn = For the general 

complete bipartite graph Kp^n+i-p see Example 3.3 below. 
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Remark 2.2. Although the explicit form of the matrix ^ = ||^fcj||i in (2-2) 
is immaterial, it is worth noting that in the unweighted case ^kn = 

1 < k < n, and one can take 

(2.8) ^ = en:=e^, 

yn 

which is known as the Fourier matrix. Clearly, there are lots of options for 
^ to be a real orthogonal matrix (rotation in 'SF with appropriate Euler’s 
angles, orthogonal polynomials etc.). 


Example 2.3. “A multiple star”. 

Consider an unweighted star-like graph Sn,p with n rays, n > 2, each of 
which contains p -|- 1 vertices, p >2. The vertices are numbered as 

{l,n-hl,...,(p-l)n-hl}, {2,n + 2,...,{p-l)n + 2}, ... {n,2n,... ,pn}, 


and the root is pn + 1, so Sn,i = Sn. Let L = Sn,p + Poo- 
As above, we put ej := ej, j = pn + 1,..., and 


(2.9) 




. ^ ^(i-l)n+q, 

q=l 


I = 


We have A(r) ej = ej-i -|- ej+i, j = pn + 2 ..., 

A(r)epfj-|-1 — 's/^^pn T Cpn-|-2) 

A(r) epn — Gpn—1 T \/U'Cpn-|-l) 

^(L) ^p(n—l)+i ^p(n—l)+i—l T ^p(n—l)+i+lj 

^(L) ®p(n—1)-|-1 ^p(n—l)+2: 


1,2, ...,p. 


i = 2,...,p-l, 


so the Jacobi subspace and Jacobi component of T are 


( 2 . 10 ) 


span{ej , 

J(r) = J({0},{a,}), a, = | 


To find the finite-dimensional component note that, by the construction, 
^p{n-i)+k £ span{e(fc_i)„_i_i,... ,efc„}. As in the above example, we supple¬ 
ment each ep(^n-i)+i to the basis in this subspace by means of the Fourier 
matrix (2.8) 


n n 

fj • ^ ^j iqj^{k—l)n+qi fn ^ ^ ^ l)n-|-g ®p(n— 

9=1 9=1 


for 1 < j < n — 1, 1 < k < p. As in (2.4) we have 


( 2 . 11 ) 


^(L)/]'^ = /f \ A(r)/f = /j') + /f ,..., 
yi(r)/f= /f+ /f^, ^(r)/f^ = /f, 


Relations (2.11) mean that the subspace FLj := span{/j^\ ..., is A(r)- 
invariant, and A{T)\7ij = A(Pp). There are exactly n — 1 such subspaces for 
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j = — 1. Finally, we come to the following canonical form for the 

adjacency matrix 

n—1 

(2.12) A(r):^ (0A(Pp))0J(r). 

i=i 


The Jost polynomial is computed in (0.11) 


^ \ ( 1 ^ 2 2 :^^ - 1 {n 

—y/nu(z) = in — l)z —---1 = — 

2 :^ — 1 

It is easy to see that the polynomial q{x) = (n- 
two positive roots 0 < xi{p,n) < Xi{p,n) = 
spectral meaning. Hence 


- 1 ) 2 :^^+^ — nz'^ + 1 
2 ^ — 1 

l) 2 .p+i exactly 

1, so the first one has the 


(2.13) crd(r) = |2cos^^|^_^|j|±(^Vxi(p,n) 


VMp, 


n 


:)}' 


SO we have p eigenvalues on the absolutely continuous spectrum [—2, 2] of 
multiplicity n — 1 and two simple eigenvalues off [—2,2]. 

Note that the weighted multiple star can be treated in exactly the same 
fashion. 

The spectrum of the multiple star was studied in [13], but no explicit 
formulae were provided. 


The problem becomes harder (in the sense of computation) if the original 
hnite star-like graph is nonsymmetric (the rays are different). 

Example 2.4. “A sword”. 

,2 

,,3 

I _^ 6 7 8 9 

11 5 Ti^2,2,oo 

,,4 


This graph (also known as T(l,2,2,oo)) can be viewed as a coupling 
T = T(l, 2,2) -|- Poo- We put = Ck, k > 6, so 


Next, let 


so that 


^(T) Cfc — Cfc-i -|- Cfc+i, 


ei + es + e^ 

65 :—- 7 =-, 64 


%/3 


A{T)eQ = Vs 65 + 67 , A(r)e5 = 


k>7. 

62 + 64 

~vr 
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We want to determine the vector 63 from the equation 

A{T)ei = 063 + 

which gives 


2 ^ ^ \/ 2 , 

- 65 = 063 + — (ei + 63 + 65), 


1 


a = 


%/3’ 


es := 


-261 +63 + 65 

V6 


and so 


^(r) 63 — ( — 266 + 66 + 6 — 2 + 66 + 6 — 4 ) — —--^=— — . 

Clearly, 63 is orthogonal to ej for j > 4 , and ||e3|| = 1 . We come thereby to 
the Jacobi subspace and Jacobi component of F 

( 2 . 14 ) /Ij = span{ej}j>3, J(r) = \/ 3 ,1,1.. . 

To compute the finite-dimensional component, we put 


62 := 


63 - 65 


61 := 


62 — 64 


^ ^/2 ’ 

so 2 l(r)e 2 = 61 , A(r)ei = 62 , and {ej}j>i turms out to be the canonical 
basis in for 2 l(r). The adjacency operator is unitarily equivalent to 


(2.15) A(r):^F(r)0j(r), F(r) = ^(P 2 ) = 

The Jost polynomial can be computed directly from relations (0.12) 


0 1 
1 0 


u{z) = ^-{1-z -3z -2z ). 

It is not hard to see that the cubic polynomial p{x) = 2x^ + 3x^ + x — 1 has 
the only real root xi, 0 < xi < 1. Hence 

1 


(2.16) 


^d{r) = { 


= \ ±1,± \/^ + 




)}■ 


3. Couplings of graphs with cycles and their spectra 

Example 3.1. “A kite”. 

2 1 



Consider the coupling of m-cycle and the inhnite path T = Cm + Poo- The 
construction depends heavily on the parity of m. 
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e,- := 


, j = n + l,...,2n-l, 




Case 1. Let m = 2n, n >2. We put as usual Ck := e^, k > 2n, and 

_ “ 1 “ ^2n—j 

V2 

Then 

^(r) e2n = \/2 e2n-l + e2n+l, 


P 4- ^n+j — ^n+j-1 + ^n+j+l 

j e2n-l — e2n-2 + V 2e2n, 

for 2 < j < n — 2, and 

44(r) — \/2 Cn + 6}^-|_2, 

The Jacobi subspace and component are 

1, 


( 3 . 1 ) Cj = span{ej}j>n, J(r) ;= J({0},{ai}), Oj = | 


\/2, j = l, 


n. 


Next, put 


e,- := 


^^7 62 n—7 


J = 1,... ,n - 1. 


' • x/2 

Then {%}j>i is the orthonormal basis in and 

y4(r)en-i = en- 2 , ^(T) ej = e,_i+ ej+i, j = 2,...,n-2, A(r)ei=e 2 , 

so finally 

( 3 . 2 ) yl(r):^F(r) 0 J(r), F(r) := ^(p„_i). 

The dost polynomial is computed in (0.21) 

—u(z) = + 2z^ — 1, 

so, by the Descarte’s rule, it has two real roots — 1 < X 2 {n) < 0 < xi(n) < 1, 
xi(n) = —X 2 (n). Hence 

(3.3) (TdiC2n + Poo) = |2cos |j|±(^xi(n) + ^ |. 

Case 2. Let m = 2n + 1, n > 1. With := e/^, A: > 2n + 1, we put 


^n-\-j • — 


%/2 


j = 1,2,... ,n. 


so 


T(r) ej = Cj-i + Cj+i, j = 2n + 2, 2n + 3,... 

7l(r) e2n+l = V2 e2n + e2n+2, ^(b) e2n = e2n-l + V 2 62^+1, 

A(r) ej = Sj-i + ej+i, j = n + 2,... ,2n - I, 

^(r) Cn+l — C^+1 T 6^+2 

(with the obvious modification for n = 1,2). The Jacobi subspace and 
Jacobi components are now Cj = span{ej}j>n+i, 

1, j + n; 


(3.4) 


•^(b) := J'({l,0,0,...},{a,}), aj = 


V2, j = 


n. 


The finite-dimensional component arises from the complement to the basis 
in the entire space 


V2 


j = 1,2,... ,n 
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and relations 

y4(r)ei =- 61 + 62 , A{T)ej = ej-i+ej+i, j = 2,... ,n-l, A{T)en = en-i, 
so finally 

(3.5) ^(r) F(r) 0 J(r), f{t) ■.= o,..., o}, {i}) 

is a Jacobi matrix of order n. 

The dost polynomial is given in (0.19) 

- V 2 {z + 1)11(2:) = 2^”+^ + 2 z ^ - 1 . 

The Descarte’s rule applied to the polynomial p 2 n+i{x) = + 2x^ — 1 

shows that there are two real roots such that — 1 < X 2 (n) < 0 < xi(n) < 1. 
Both of them contribute to the discrete spectrum of the Jacobi component 

o-d(J(r)) = |xi(n) H-^ ,X2(ra) H-. 

t xi(n) X2(n)J 

To find the spectrum of F{T) in (3.5) we expand the characteristic deter¬ 
minant over the first row 

det |F(r) + x\ = det |^(Pn) + a:| - det |^(Pn-i + x\ = Un - Un-i , 
where Un is the Chebyshev polynomial of the second kind, 

, , sin(n + l)t , rr ^ \ ^ 

Un[COSt) := -—- , Un[COSt) - Un-l[COSt) = - 7 - . 

sm t cos I 

So the spectrum of F{T) is 


r(F(r)) = | 2 cos^^r . 

^ ^ I 211 + 1/7=1 


Actually, we can combine the latter two results in a unique way. 
Proposition 3.2. The spectrum of the graph T = Cm + Too fs 
Cj(Cm + Too) = [~2, 2] U (Trf(Cm + Tqo) 
with the discrete component 

™ . r, 27r/c'i' 11 ... ... . \m 


(3.7) cJd(Cm + Too) = |2cos-^|^_^[J{Ai(m),A2(m)}, ^ — , 

where 

Xj ( m ) = Xj ( m )-\ - 7 — r , 7 = 1,2, — 1 < X2(m) < 0 < xi(m) < 1 

Xj [mj 

are real roots of the polynomial p{x) = x™' + 2x^ — 1. 

Example 3.3. “The complete bipartite graph”. 

Let G = Kp^n+i-p be the complete bipartite graph of order n + 1 [7, p.l6], 
T = iLp,„+i_p + Too. Define 

Si := 6i + . . . + 6p, 5*2 := 6p+i + ... + Cn, 

and put 


67 := 6 j, j = n + l,n + 2 ,..., 6 n := 


^n—1 • — 


'U — p 
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Then A(r) Cj = ej-i + ej+i for j > n + 2 and 

^(r) Cn+l — <51 + (in+2 — '\/p(^n T (^n+2- 

For the bipartite graph we have A(r) = 52 + Cn+i for 1 < /c < p, and so 

yl(r)en = ^p( 52 + en+i) = ^Jp{n-p)en-l + y/pen+i- 

Vp 

Next, since yl(r)efc = 5i for p + 1 < A: < n, then A{T)en-i = \/p{n — p) e^, 
and the Jacobi component is 

(3.8) £,/= span{ej}j>n-i, J(T) := J{{0},{y^p{n - p),y/p,lA, ■ ■ ■})■ 

To find the finite-dimensional component, note that G spanjcfc}^^^, 
Cn-i G span{efc}^^p_,_;^, so we amplify these two vectors to the orthonormal 
base in the corresponding subspaces by means of the Fourier matrices ^ = J^p 
and T] = Fn-p 


p 

^ y ^kj^k-i J ~ 1) • • • 1) 

k=l 
n 

dj := ^ ^ j = p, ■ ■ ■ ,Tl 2, 1 

k=p+l 

We have ^(F) = 0 for 1 < /c < n — 2, so 


1 




'ypti 

1 

y/n-p 


Y1 


k=p+l 


(3.9) A(F) ~ F(F) 0j(r), F(r):=o„_.. 

The Jost polynomial is given in (0.16) 

—p^/n — pu{z) = {p — l)z^ + (p(n — p) + p — 2)z^ — 1 = q{z‘^)- 

It is easy to see that the quadratic polynomial q has the only root xi(p,n) 
in (—1,1), 0 < xi{p, n) < 1, so 


(3.10) cj(i(iiip^„+i_p J-Poo) 


{oi”-2>} y {± (vAFJN) + }. 

^ Vxi(p,n)^J 


For p = 1 we have Example 2.1. 

Example 3.4. “A propeller with equal blades”. 
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The graph T is composed of two equal cycles C 2 n+i having one common 
vertex, with the infinite path attached to it. We put eu = ek, k > An + 1, 

- + ^ 2 n-i + e 2 n+i+l + e^n-j , _ n i „ i 

^4n—j •— 2 ’ ^ 

Then 

^(T) G4n+1 — 2e4}^ T B4n+2^ 2 l(r) G4n — 64 ^ 1—1 T 2e4n-|_i, 

y4(r) G4fi—j — G4ji—j — l T 64 n_j-|_l, j — 1, 2, . . . , U 2, 

2 l(r) 63^+1 = esn+i + e3„+2- 


So the Jacobi subspace and Jacobi components are now Cj = span{ej}j> 3 n+i, 

1, j + n; 


(3.11) 


J(r) := J({l,0,0,...},{aj}), aj = 


2, j = n 


(cf. (3.4)). 

The finite-dimensional component arises by putting 


fr-= 

93 ■= 


hj : = 


V2 

^3n—j + l ^Sn+j 

71 ’ 

®n—j+l T ^n+j 63^_j_|_4 G^n+j 


3 = 1, 


,n. 


Each of them satisfies 


Vi = - 91 + 92 , A{T)yk = yk-i+9k+i,k = 2,...,n-l, A{T)yn = yn-i, 

so each of linear spans is n-dimensional and ^(T)-invariant subspace, and 
(cf. (3.5)) 

i"(r) = 0^({-i,o,...,o},{i}). 

k=l 

The rest is the same as in Example 3.1. The Jost polynomial is 
-2{z + l)u{z) = 3^2”+^ -h Az^ - 1, 

and it has two real roots such that —1 < X 2 (n) < 0 < xi(n) < 1. Both of 
them contribute to the discrete spectrum of the Jacobi component 

crd{J{+)) = {xi{n) H-^ ,X2(n) H-. 

1 xi(nj X2(ra)J 

The spectrum of E(r) is given in (3.6), but now each eigenvalue has multi¬ 
plicity 3. 

For spectral properties of finite propeller graphs see [15]. 


Remark 3.5. If the blades are C 2 n) the argument goes through in exactly 
the same way (see Example 3.1, case 1). Moreover, the algorithm applies 
equally well to propellers with p equal blades for p > 2 with the infinite path 
attached to their common vertex (“flowers”). 
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The first step is standard 


Cfc := efc, A: > 5, 64 := 




A{T)ek = ek-i + ek+i, k> 6 , yl(r) 65 = 2e4 + ee, 
yl(r)e4 =- - -h 265. 

We find 63 from 063 = 4 l(r)e 4 — 6 E 4 — 265 , where the parameters a,b are 
determined from the orthogonality ( 63 , 64 ) = 0 , and normalization ||e 3 || = 1 . 
It is easy to see that 

«= 2 ’ ^= 2 ’ ^- 2 -’ = 2^3 + 2^4 + 265 . 

Hence 

Af-r\'- _ei+e4_ 1 ^ 1 ^ 

^(r) 63 —- - - “ “2 2 

so the Jacobi subspace and Jacobi components are 

(3.12) £j = span{e,}j> 3 , J(r) := ^ , 0,.. .|, , 2,1,.. .|^. 

Next, put 

61-64 _ 62-63 _, 4 

ei := ^ , 62:= ^ ^>H(r)6i=62, yl(r) 62 = 61 - 62, 

and so the canonical from is 


(3.13) 


H(r):^F(r)0j(r), F(r) := 


0 1 


The Jost polynomial is now 

—u{z) = + 2>z^ + 3z^ + 2 ; — 1, 

and it has two real roots —l<X2<0<xi<l. Hence, the discrete 
spectrum is 


^d(r) = { 


-1± V5 1 1 1 

--- ,xi H-,X2 H- >■ 

2 Xi X2^ 
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Example 3.7. “A wheel”. 

The wheel Wn is a graph consisting of a cycle on n vertices, 1, 2,..., n, 
n > 3, and the vertex n + 1, adjacent to each of 1, 2,... , n (cf. [1, p. 49]). 
Consider the coupling T = Wn + Poo with the path attached to the root 
n + 1. We put 

ej:=ej, j>n + l, 

and so 

A(r) ej = ej-i + Cj+i, j >n + 2, 

yl(r) Bn — ‘2‘Gn T P ^ 65 , 

SO the Jacobi component is 

(3.14) Cj = span{ej}j>n, J(r) := J({2, 0,0,...}, {\/n, 1,1,...}). 

To find the finite-dimensional component, consider the adjacency operator 
A{Cn) acting in C"^. Let 

{/l,/2,...,/n}, /,GC-, /„ = ^{1,1,...,1} 

Vn 

be the orthonormal basis of its eigenvectors. Take ej := {fj,0, 0,...}. It is 
a matter of a simple computation to verify that 

= {A{Cn) fj,0,0, ...} = 2 cos ^^^^ ej, j = 1 ,... ,n - 1 , 
so the canonical form is 

(3.15) A(r) ~ F(r)© J(r), F(r) := diag|2 cos — ,..., 2 cos ~ 

In n J 

The dost polynomial is computed in (0.15) 

— \/nu{z) = (n — l)z^ +2z — 1. 


Cn • — 


n 


k=l 


A(r) Gn+l — \/n Gn P Gn +21 


It has two roots 

X2{n) = 


\/raP 1 

n — 1 


< 0 < xi{n) 


y/n-l 
n — 1 ’ 


such that xi{n) G (0,1) for all n > 3, and for X 2 (n) we have 
X 2 {n) < —1, n = 3,4; X 2 {n) > —1, n > 5. 

Hence, the discrete spectrum is either 

{ OikiT ^^ I I r 1 'i 

or 

o-d( Wn PPoo) = (2 cos I M|a:i(n)H-^,X 2 (n)H-n > 5. 

L n ) k=i ^ t xi(n) X 2 yn)) 
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4. Graphs with several tails 

Given a finite graph G, we can attach several infinite paths to one or 
several vertices of G. Although such graphs are not exactly couplings in 
the sense of Definition 0.1, the algorithm suggested in Theorem 0.6 works 
in this case as well (see Remark 0.8). 

Example 4.1. “A cycle with a double tail”. 

Given the cycle C 2 n+i with vertices 1,2,... , 2n +1, we attach two infinite 
paths to the last vertex. Denote this graph by T = C 2 n+i + 2Poo. 



We proceed as in Example 3.1, case 2, (see Remark 0.8) and put 
en+j ■= -^^ ej := -^^j = l,2,...,n. 


and 


\/2 ’ V2 

C2n+fc T C2n+fc+l 


^2n-\-k • — 


V2 


, A: = l,2,.... 


Clearly, the adjacency operator on the A(r) -invariant subspace C = span{ej}j>i 
acts exactly as the one in Example 3.1, so the spectrum of A(r) \C is known. 

The rest is obvious. The orthogonal complement Q L is spanned by the 
system of vectors 

T ^2n-\-k ^2n-\-k-\-l 7 

Ik :=- ^ -, A; = 1,2,, 

and A(r) = Jq. So 

(4.1) ^(r):^F(r)0j(r)0Jo, 

F{T), J(r) are given in (3.5). 

Example 4.2. “A cycle with two tails”. 
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Given the cycle C 2 n+i with vertices 1,2,... , 2n +1, we attach two infinite 
paths to the last two vertices. Denote this graph by F = C 2 n+i +lPoo + lPoo- 
We put Cl = en, 


_ en+j-l + Cn-j+l _ o o „ 

, J — Z,, O, . . . , / 6 , 

~ _ e2j-2 + e2j-i 




, j>n + l. 


The subspace L ;= span{ej}j>i is yl(r)-invariant and 

y4(r)ei = V 2 e 2 , 4 l(r) 62 = ^2 ei + 63 , A{T)ek = ek-i + ek+i, 3<k<n 

^(T) Cn+l = Cn + e.n+1 + 6^+2; ^(r) &k = ^fc-l + ^k+li k > n + 2 


So ^(r)| C = J{{bj}, {%}) with 


(4.2) 




1 , j = n + l, 
0 , j ^n + l. 


Next, put 


a,- = 


V2, j = l, 
1 , j / 1 - 


^ Cn+j Cn-j 


1 < i < — 1, /ij : = 


_ e2i - e2j+i 


j > n. 


^ ■ V2 

It is clear that Ti = span{/;,j}j>i is ^(r)-invariant, A{T)\7i = Jq, and 
C^BT-L = ^. Finally we have the following canonical form 

(4.3) ^(F):^ J({6,},{a,})0Jo, 


J{{bj},{aj}) is in (4.2). 

The dost polynomial is computed in (0.20) 

— \/2u{z) = + z — 1, 


and, by the Descarte’s rule, we see that u has one real root 0 < xi{n) < 1. 
Hence, the discrete spectrum is 

o-d(r) = xi{n) H-, 

xi(n) 

and the absolutely continuous spectrum has multiplicity 2 . 


Surprisingly enough, the case when p > 1 infinite rays are attached to 
each vertex of a finite graph G with vertices {1,2 ,... ,n} (so there are pn 
such rays altogether) is easy to handle, and the spectrum of such graph can 
be found explicitly in terms of the spectrum of G. Denote such graph by 
F = G + Poo (p) ■ We label the vertices along the rays (off G) as 

Ri = {ruj + i}]Li, rrij := n + {j - l)pn, i = 1,2,... ,pn. 

Theorem 4.3. Given a finite graph G with vertices 1, 2,..., n and the spec¬ 
trum cr{G) = let r = G + Poo(p), p G N. Then the adjacency 

operator H(F) is unitarily equivalent to the orthogonal sum 

n 

A(r) = 0j(A„Vji)0( 0 j„), 

j=l i=l 

J{>^j,Vp) ■= «/({Aj,0,0,...},{ Vp, 1,1, •••})• 


(4.4) 
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The spectrum of F is 

n 

(4.5) a{r) = |-2,2]UU 

i=i 

Proof. Put 

1 ^ 

gj{k) :=— '^eraj+(k-i)p+i, k = l,2,...,n, jeN, 

1=1 

so {gj{k)}j^k is the orthonormal sequence in It is easy to see from the 
way of labeling of the vertices that the adjacency operator 44(r) acts on 
these vectors as 

M'^)gjik) = gj-iik) + gj+iik), j = 2,3,... 

A{'^)gi{k) = Vpek + g2{k). 

Next, by the definition 

(4.7) A{T) ek = {A{G)e^^\o,0,...} + ^gi{k), k = l,2,...,n, 

where is the standard basis in 

To construct the canonical basis we invoke the eigenvectors of 44(G) 

n 

ft 'c = iic,iii;,tei 

q=l 

is a unitary matrix, so A{G)fjf^^ = \kf^\ k = l,...,n. With such a 
notation at hand we put 

n 

eo{k) := ^CqkSq = {fj!"\o,0 ...), 

<?=! 

n 

^jik) ■.= '^<fqkgj{q), j G N, /c = l,...,n. 

<?=! 

It is clear from (4.6) that 44(r) ej{k) = ej-i{k) + ej+i{k), j > 2, 

n 

A{T)ei{k) = iVP(^q + 92iq)) = Vpeoik) + e 2 {k), 

<7=1 

and 

n n 

44(r) Fo(fc) = Cqk ^(r) e, = (J] Cqk A(G)e(-), 0,0,.. .) + ^Fi(fc) 

q=l q=l 

= Afceo(A:) + ^ei{k). 

Hence the subspace C := span{ej(A:)}i<fc<nj=o,i... is 44(r)-invariant and 

n 

41(0 |£ 44 0 J (Afc) \/p)- 

fc=i 
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It is not hard to see that on the orthogonal complement 

(p-l)n 

A{r)\feCc^ 0 Jo. 

k=l 


The proof is complete. □ 

In the case p = 1 this result is proved in [18]. 

To find the discrete spectrum of J{^k,y/p) one has to solve the dost 
equation (0.15) 

(4.8) {p - l)x^ + Xkx -1 = 0, 

pick its spectral roots, that is, the roots in (—1,1), and then take their 
Zhukovsky images. The simplest case is p = 1, when each with |Afc| > 1 
generates one eigenvalue A^ + A^^ of A(r). For p > 3 an elementary analysis 
of equation (4.8) shows that there are two spectral roots for |Afc| < , 

and there is one such root otherwise. 

Example 4.4. Consider the case p = 2. Equation (4.8) has the roots 

Afc ± yjxl + 4 

x± = - - -, x_ < 0 < x+. 

It is easy to see that there is exactly one spectral root for Xk / 0, and no 
such roots for Xk = 0. The discrete spectrum now looks as follows 

o'd(r) = IJ signAfc ^Jxl + A. 

Afc JO 

We complete the section with two examples wherein two-sided Jacobi 
matrices arise naturally. 

Example 4.5. “A double infinite star graph”. 

Denote by Sk,oo the infinite star-like graph with k > 2 infinite rays em¬ 
anating from the common root (cf. Example 2.3). The main object under 
consideration is the coupling of two such graphs with the bridge connecting 
their roots, so T = Sp^oo + Sq^oo- 

We label the vertices of the “right” star Sp^oo along the rays by positive 
integers as 

{1, 2,p-|-2,..{1, 3,p-|-3,, {l,p-|-1,2p-|-1,...}, 

so the root is 1. Similarly, we number the vertices of the “left” graph Sq^oo 
by negative integers, so the root is —1 and the vertices along the rays are 

{-l,-2,-q-2,...}, {-l,-3,-g-3,...},..., {-1,-q - 1,-2q - 1,.. 

The underlying 1“^ space is £‘^{Zq), Zq := Z\{0}, and the standard basis is 
{cjljeZo- 
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We construct an orthonormal system by 


Cl ei, 

Cfc := 

— ^e(fc_2)p+i, k>2; 

V i =2 

1 

J<-2- 

e_i := e_i. 

+j := 

is easy to see that H(r) 
y4(r)ei = e-i + ^/pe 2 , 

acts as 

yl(r)e^_i = ^e^_2 + 61, 

H(r) 62 = y/pei +63, 


H(r) e^_2 = e_3 + y/qe-i, 

2l(r) Cfc = Ck-l + Cfc+I, 

k > 3, 

A{T)e-j = e-j-i + e-j+i 


Hence the subspace L = span{ej}jgZo is H(r)-invariant and 

{ Vp, J = 1; 
i = -1; 

1, 3 + ±1. 

To supplement the system {ejjj^Zo to the orthonormal basis we proceed 
in a standard way 


^0 


p+i 


^ ^ ^ t) • • • ^ 1)2,.. 


i=2 

9+1 


i=2 


^ ^ Vi—1,1 ^{j+l)q—ii ^ 1) • • • ) 9) J 1) 2, 

Finally, 

(4.10) 


IlCrs II — Tp; 


?rs|| — J'q- 


p+q-2 


H(r) J(p,g)0( 0 Jo 


i=l 

To find the discrete spectrum we apply Example 0.5 with a_i = ^/q, 
= v^) so characteristic equation (0.35) looks 

Qiy) = (p - i)(^ - i)y^ - (p + g - l)y + 1 = 0. 

It is easy to see that this equation always has tow positive roots 0 < y- < y+. 
In the cases min(p, q) = 2 or p = q = 3 we have > 1, so only the first 
root has the spectral meaning. Otherwise, 0 < y_ < y+ < 1, so both of 
them are the spectral roots. Hence, 


or 


(^diJip,q)) = ±(^/y^ +mm{p,q)=2 or p = q = 3, 

= {±(vw+ ^).±(vw+ ^)} 


for the rest of the values p,q >2. Here 

p + q-l± \/(p^^gP”+2(pT”gy^^ 


y± = 


2ip-l)iq-l) 
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Note that if the bridge has weight d, the above argument ends up with 
the Jacobi matrix 

{ y/P, j = 1; 

y/q, J =-1; aj = 1, 

d, j = 0, 

for which the spectrum can be found explicitly (see Introduction). 

The example was studied in [18] by using a different model based on 
one-sided Jacobi matrices (0.5), with no explicit expressions provided. 


Example 4.6. “Infinite regular trees”. 

The Bethe-Cayley trees of degree d G N (or Bethe lattices in physical 
literature) are among the most notable infinite trees. For the introduction 
to the spectral theory on such trees see [21, Chapter 10]. We consider here 
the coupling of two copies of B^ with the bridge connecting their roots, 
r = Bd J- Brf. Note that T is an infinite regular tree of degree d + 1. We 
number the vertices of T as in the previous example (the vertices of the 
“right” copy by positive integers, and of the “left” copy by negative integers). 
So, the vertices of fc’s generation in the right copy (the root forms the 1st 
generation) are labeled with 


{h-i + 1 ) h-i + 2 ,..., Ik}, Ik ■= d^ = ——Y ) k gN. 

i=i 

Again, the underlying Hilbert space is ^^(Zq) with the basis {ejljg^o- 
Put e±i := e±i. 


Gk 


1 

Vd^~^ 


£ ei, 

i —1 H“ 1 


^—k • — 


1 

'/dF^ 



J—h-1+1 


k = 2,3,.... 


The operator A(r) acts as 

A(r) ei = e_i J-\/de 2 , A(r) e_i = \/de _2 J-ei, 
A(r) ek = Vdek-i + Vdek+i, \k\ > 2 . 


So the subspace £ = span{ej}jgZo is A(r) -invariant and 

(4.11) A(r)|£ = J(^) = ^/dJ({0},{a,},6z), = 

The structure of operator A(r) on the orthogonal complement is the same 
as in the case of B^ (see [3], [21, Theorem 10.2.2]), so finally 

(4.12) A(r) J(r) 0 (0 ^/d Jo) , = I j ^ f 

i=i I ’ ■ 

To find the discrete spectrum of A(r) note that J(r) (4.11) is the subject 
of Example 0.4. It follows that the discrete spectrum of J(r) is empty, and 
so 

(4.13) u(r) = [-2\/d, 2\/d] 

and is pure absolutely continuous of infinite multiplicity. 
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Remark 4.7. Since Bi = Pqo) the question arises naturally whether it is 

possible to find spectra of the couplings F = G + The above algorithm 

works equally well in this situation and leads to the canonical form 

A{G + B(i) ~ A{G + Poo) ©(© VdJoj . 

i=l 
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